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Abstract 

The paper gives a screw axiomatics of rational mechanics, namely: 

1. introduces the main measures of mechanics: the mass measure, the scalar and (vector) screw 
measures of motion, the (vector) screw measure of impressed action, the increment velocity 
of the vector measure of motion, the (vector) screw measure of constraint action; 

2. postulates the (stronger) local integral form of conservation law for the vector measuare of 
motion (fundamental principle of dynamics), and 

3. defines the central concept of rational mechanics - mechanical system being realized in the 
form of all classical mechanical systems ( mass-points, rigid bodies, continua, point-bodies, 
etc.). 

The presentation is based on new notions of vector calculus - homogeneous and inhomogeneous 
vector and tensor slider-functions and screw measures. 

PACS numbers: 45.20.D-, 46, 83.10.Ff, 47.10.ab, 83.10.Gr. 

Keywords: classical mechanics, continuum mechanics, constitutive equations, measures, foundations of me- 
chanics, screw theory. 



1. Elements of rational mechanics 

Due to [1] 'Rational Mechanics is the part of mathematics that provides and develops log- 
ical models for the enforced changes of place and shape we see everyday things suffer. 
. . .Mechanics does not study natural things directly. Instead, it considers bodies, which are 
mathematical concepts designed to abstract some common features of many natural things. 
One such feature is the mass assigned to each body. Always, a natural body is at any one 
instant found to occupy a set of places; that set is the shape of that body at that instant. . . 
The change of shape undergone by a body from one instant to another is called the motion 
of that body. . . .motions of bodies are conceived as resulting from or at least being invari- 
ably accompanied by the action of forces. . . Mechanics relates the motions of bodies to the 
masses assigned to them and the forces that act on them. Bodies are encountered only in 
their shapes. Masses and forces, therefore, can be correlated with experience in nature only 
when they are assigned to the shapes of bodies '. 

It is easy to see that the statements quoted above lead to the contradiction as the working 
tools of mechanics are connected with the measures carried over to body shapes, instead of 
bodies. This contradiction is assumed to be removed by means of carrying over the forces 
acting on bodies to those shapes 'in some special way' PQ. 

Furthermore [ilj: . . . the basic categories of mechanics have dual character: on the one hand, 
in essence such concepts as 'force' and 'mass' are abstract mathematical categories. On the 
other hand, they appeal to natural objects with which human practice deals: 'force' represents 
a measure of object interactions, measured by any physical means; 'mass' - quantity of the 
substance containing in an object'. 

In this connection, it is of importance to notice that, due to this dual character, 'body' is 
considered as a natural thing or as the mathematical concept 'designed to abstract some 
common features of many natural things', but not only - mass. 

It is the concept of body (mechanical system) that we shall refine as a mathematical category 
in this paper. 

Below, we shall use the set R of all real numbers and n-dimensional affine space A n modeled 
on n-dimensional vector space V n . 
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1.1. Slider— functions and screw measures 

We shall introduce the notions of vector and tensor sliders being a key moment for the 
theory of mechanical systems which is represented below. 

Due to the Great Soviet Encyclopedia, v. 5 (Moscow: Soviet Encyclopedia, 1971), 'screw 
calculus is the section of vector calculus in which operations over screws are studied. Here the 
screw is called the pair of vectors {p, q}, bounded at a point O and satisfying to conditions: 
at transition to a new point O' the vector p does not change, and the vector q is replaced 
with a vector q' = q — 00' x p where x means cross-product. The notion of the screw is 
used in the mechanics (the resultant / of a force system {fl} and its main moment m form 
the screw {/,m}), and also in geometry (in the theory of ruled surfaces)'. 

Note that we do not support the idea to use the name 'torser' from the French 'torseur' 
instead of 'screw' [2]. 

Let us add more details to this notion. 

1.1.1. Slider vector— functions. Assume that there are vectors p x and q x bounded at a 
given point x G A3, and at any point y G A3 

Py=Px, q y = q x + ry X p x (1) 
where Vy X is the spin-tensor generated by the vector f yx = yi. 

Definition 1. The field l Px > qx = {p yi q y ,Wy G A 3 } is called slider vector-function or, briefly, 

def _, 

slider while l Px ' Qx = {p y ,q y } is known as reduction of the slider w.r.t. the reduction point 
y e A 3 . 

A slider is called homogeneous if q x = 0. In this case we shall use the notation l Px . 

If one marks coordinate columns of vectors in a Cartesian frame £ = (Oq, e ) (Go is its origin 
and e is its base) with the superscript °, then l Px ' 9x '°={Py,qy,Vy G A 3 } is the coordinate 
representation of the slider l Px ^ x . In order to apply the matrix tools one may use the 
following coordinate columns lP x >i x ' wr >° = col{p°,g°} and lP^i^ tw '° = col{g°,p°} known as 
wrench and twist (in the space E 6 ), respectively. 

1.1.2. Screw measures. Let 0-3 be a-algebra of subsets in A3. Introduce the following 
Borel measure 

fi(A) = Hac(A) + Hpp(A), Aea 3 
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where // ac (A) is the absolutely continuous component w.r.t. Lebesgue measure /X3 and (i pp (A) 
is the pure point (discrete) component presented as n PP (A) = Ylk^k f° r points x k e A whose 
are called pure, the others being called continuous [3]. 

Definition 2. Let \ A be the characteristic function of A. The Lebesgue-Stieltjes integral 

ir(A) = J X J Px ' q ^(dx), A G a 3 (2) 

is called signed screw measure or screw. Signed measure is the generalization of the notion 
of measure by allowing it to have negative values 

We shall use this name for surface Lebesgue-Stieltjes integrals, too. 

The screw measure is a screw in the sense of the Encyclopedia definition (if {p, q] d = 7Tq(A) 



def 

at the point O then there is n'o(A) = {p, q'} at the point O' where q' = q — 00' x p). 

Screws generated by homogeneous (inhomogeneous) sliders will be called homogeneous (in- 
homogeneous). 

Define the triple of orthogonal unit vectors e±, e<i and e*3 in the 3-dimensional space V3. Let 
us introduce 6 screws such that at a point y G A 3 their reductions are defined as follows 

-0— (*)•-(*)•**-©•-(*)•- (i 

where o G V3 is the null vector. 

As any screw is defined in the unique way by its reduction at some point, these 6 screws 
generate a basis of 6-dimensional vector space (see also [2]). 

1.1.3. Multiplicative groups of motions. In A3 let us introduce the Cartesian frames 
£ p = £{O p , e p ) and £k = £(Ok, e&) with the origins O p and O k and the bases e p and e^ where 
p and k are naturals. 

Define the rotation matrices Co, p and C p ^ defining orientations of the Cartesian frames £ p 
and £k w.r.t £ and £ p , respectively. Then C 0tP C p ^ = Cq^ and for any free vector A there 
are the following relations 

A = Co iP X p , X p = C Pi kX k 

where A , A p and A fc G E3 are the coordinate columns of the vector A in the bases eo, e p and 
e k . Hence we have also A = C 0) fcA fc . 



Introduce the radius-vectors f x and f p>x G V 3 of a point x G A 3 w.r.t. the origins Oq and O p , 
respectively. Define d 0)P = r x — f P)X . Then we may represent the relation f x = d 0)P +f PiX G V 3 
in the frame £ as r x = d[j + C 0>p r^ x . Let the point x be immovable in £ p then r v px is time- 
constant. 



With differentiating the last relation we have v x = Vq^+Cq p r^ where v x = r x is the velocity 
of x and vo, P = d^ p is translation velocity of £ p w.r.t. Oq, respectively. Here to honor Newton, 
we use the superscript " for full derivatives by t, e.g., for any function A = A(x(t),t) we have 
A' = §- t A + (divA)x-(t). 



Hence 



For any vector / 



h 



' j o,p + C Pi oC 0tP rp >x 



(3) 



G E 3 introduce the cross product matrix 



r 



dcf 



o -h h 

h -f x 
-/ 2 /i 



(4) 



Thus we may define (in £ p ) [IT] : 



the coordinate column d^ p G E 3 of the translation vector of £ p w.r.t. £o in £p! 
the coordinate column Vq G E 3 being known as quasi-velocity of the translation 



of £„ w.r.t. £n in £ 



the cross product matrix oJq* 



dcf 



C Pj oCo » where the triple wj GE 3 is known as 



quasi-velocity of rotation of £ p w.r.t. £q (calculated in E i 



pj- 



The algebraic quantity u p p corresponds to the geometrical one, i.e., the instantaneous an- 



gular velocity u; ,p G V 3 . It defines the rotation axis of £ p . 



Introduce following matrices 



0,f 



0,p 



O 



o a 



0,p 



0,p 



d 0x 

u 0,p 



O 
I 



d px 

a 0,p 



O 
I 
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Theorem 1. [5j A given screw n, -Rq 



wr,0 



LZ P K r,p where n o 



wr,0 



and 7r™ r,p are wrenches of 7r 



computed in the frames £q and £ p , respectively, the matrix Lq^ has the representat 



ion 



jwr 7-)0 /^y® r-)P 

n 0,p ~ -^O.p^O.p — °0,p- ty 0,p 



(6) 



and belongs to the multiplicative group £ wr (TZ,6) such that 

px 



t wr. t wr /Rwr 

L Q,p — ^O,p^0,p 



lTrwr t wr rf\Wr 



px px 
V 0,p u 0,p 



ywr 
0,p 



< O 

v 0x uj 0x 
u 0,p u 0,p 



(7) 



Proof. Representation ^ follows directly from the screw definition. 

Consider the case where L^ p = C® p Dl p . Then from g follows that = C®; p D p Qp + 

\ U 0,p U 0,p + U Q,p)^Q,p U ^p — 



<® n p - — (n®- r>p dp r 1 ® _i_ r y ® n p - n p r 1 ® ^® n 

O.p-^p ~~ ^0,p- ty 0,p- L/ p,0'^p,0 "r < ^0,p- Ly 0,p- ty p,0 l ^p,0>' l - y 0,p- L/ 



P 

0,p 



iT/wr T wr 
^0,p-^0,p- 



" uir 
v p,fc 



The matrices of the kind L^ p = Cq p Dq p form a group because there are Lq^L^ 

r~i<S> /n® /-rig) Y)P /-rig) r)k _ #~r® r)k r)k _ /nr® p,k _ j wr „„J jwr-\ _ / f~t% n p \_i _ 
^^p^p,k^k,p U ^p^p,k U p,k — 0,P PM ~ °0,fc W 0,J! — n 0,k anQ -^0,p ~ l^Cp-^CpJ ~ 

/rpp \_l^-y(g),T _ rpp ^0 _ .-rig) rpO _ T wr 

\ 1 0,p) o 0,p ~~ 1 p,0 Ly p,0 ~ Ly p,0 1 p,0 ~ ^p,0- 

Consider the case where Lg£ = 7>° p C ® p . Th en from g follows that L^ p = D$ p C® p + 

nO /-r(g>. nO /-rig) //-rig) /-rig>. _i_/-rig> nO piO. /-rig \ jwr^wr 

ly O,p u 0,p ~~ ^0,p u 0,pl u p,0 u 0,p " t_L 'p,0- ty p,0- ty 0,p L/ 0,pJ ~~ -^0,p^0,p- 

The matrices of the kind = 7?g p C® p form a group because there are Lq^L^j. = 

rpO /-rig) rpp sy® _ rpQ rpQ ,-r(gj _ p.Q ,-r(g> _ T wr Q „J T Wr -1 _ j rpQ .-rig) \_1 _ .-rig: /rpO \ — 1 _ 

1 0,p^ J 0,p 1 p,k ( - J p,k ~ 1 0,p 1 p,k^'0,k ~ • L, 0,fc L/ 0,ft ~~ ^0,fe dnu -^0,p ~~ I-'O^OjjJ ~~ 0,pJ ~~ 

rpO /~t®,T ,-rig) ry-ip .nflg) r wr 

( - / p,0- t p,0 , ~ y 0,p u 0,p ~~ - t p,0 ( -'p,0 ~~ "^p.O- 

As result we have the following relations 



r uir. 

^0,p 



T wr (f) wr 
-^0,p^0,p> 



^0,p 



^O.p-^O.p 



(8) 



being matrix-functions of (JJ^, i^, and 



UJ, 



px 
0,p- 



The similar statement 7rf, w ' p = L 1 ^" n^' is true for twists where we have the matrix 



T tW 

^0,p 
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-^0,p 







/ 
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(9) 



belongs to the multiplicative group C tw (K,6) such that L l w p = % w p L^ p , % w p = -&™ p ' T and 

wr,T 



T tw. rf\tw T tw ffttw 

^0,p — ^0,p-^0,p' ^0,p 



O.p 



Note that in contrast to the groups of motions in the 3-dimensional space the groups 
£ wr (1Z, 6) and C tw (1Z,Q) are multiplicative. 



1.1.4. Slider tensor— functions. The slider notion is based on the pair of vector-functions 
p x and q x . That is why these sliders are called vector ones. If we replace these vector- 
functions with tensors V x and Q x of II rank, then the corresponding sliders will be called 
tensor ones. 
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1.2. Main concepts and structures of mechanics 

In what follows, we shall use Galilean space-time [6] introduced as the quadruple G = 
{V 4 , A 4 , t, g} where 

1. t: V 4 — > Vi is a surjective linear mapping called time one, and 

2. g = (•, ■) is an inner product on ker{r} (= V 3 ). 

The points of A4 are called world points or events. The number r(b — a) is called time 
interval between events a and b G A 4 . These events a and b G A 4 are called simultaneous if 
r{a — b) = 0. The set of simultaneous events forms 3-dimensional affine space A 3 in affine 
space A 4 . 

The inner product (-, ■) (in Galilean space-time) enables one to pass from the space V3 to 
Euclidean space E 3 with the norm \\x\\ = -J (x, x) and to introduce Cartesian frame £q in 
A3 (with the origin Oq G A3). 

The basis e = (e ,i, e 0j2 , e , 3 ), e ,i = col{l,0,0}, e , 2 = col{0, 1,0}, e , 3 = col{0,0, 1} is 
called canonical. Free vectors and their coordinate columns in the canonical basis eo are 
known with the same name vector as elements of the vector spaces V3 and E3. In what 
follows, one uses the number space E 3 as representation of V 3 . 

Any set T C R may be used for parameterization of the image of r with cr-algebra a t of 
subsets in R. Values of parameter t G T are called instants. We shall assume that there is 
defined cr-algebra a t and the Lebesgue measure fi(dt) on the set T. 

Let us introduce following notions [TJ. World-line is a curve in A 4 whose image in A 3 x T 
associates one point x(t) G A 3 to each instant t G T. A collection of non-intersectional 
world-lines forms world-tube. Here 'Intersections of world-lines represent collisions or the 
creation or destruction of bodies or elements of bodies. In specific mechanical theories such 
intersections are usually excluded (the principle of impenetrability) altogether or allowed as 
exceptional cases subject to specified conditions' pp. 

Henceforth we shall name some world-tube A C A 4 as universe. As in probability theory 
[7] , the universe is separately specified for every mechanical problem. 

A given world-tube Ac A, the world-tube A e = A \ A is called environment of A in the 
universe. 

The universe A defines the family {A t C A3, t G T}, for any world-tube Ac Awe having 
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the family {A t C At, t G T}. We shall assume that the Borel measure introduced above is 
time-invariant on the sets A t and, if a point x(t*) of any curve {x(t) G A t , t G T} is pure (or 
continuous) at some time instant £*, all points of this curve are also pure (or continuous). 

Let us introduce kinematical and dynamic structures in the universe. 



For each point x(t) G A t , t G T, the radius-vector f x (t) = (O ,x(t)) is called position of the 

def 

point, and a vector v x = v(x(t),t) = f x \t) is called its velocity w.r.t. O at instant t G T. 

We shall call mass the scalar measure A4(A t ), being continuous w.r.t. p(dx). Due to Radon- 
Nicodym theorem the measure may be represented as the following Lebesgue-Stieltjes inte- 
gral 

M(A t ) = J x At p x n(dx), A 4 C A t (10) 
where p x is the mass density. 
Let us define the following scalar measure 

£(At) = J X At k x p(dx) (11) 

with the density k x = \{v X) p x v x ) . 

Introduce the following vector and screw: 

d 



-K = PxV x , V{A t ) = J x A J Px p(dx) 

Due to j5] we shall use the notion of bi-measure: a vector-function $(•, •), which is defined on 
o",3 x cr 3 and a screw measure of the kind ^ by each of arguments, is called screw bi-measure. 

A bi-measure <P(A, B) is called skew if @(A, B) = -<P(B, A) for any A and B G a 3 . 

Let a skew screw bi-measure $(A t , A^) be homogeneous. By definition it is the screw w.r.t. 
every argument. That is why there exists the slider P x such that the bi-measure coincides 
(by the first argument) with the following screw 

J-(A t ) = I X J f *p(dx) = f <2>(A t ,A t e ) (12) 



and «P(A|,A t )=-^(A t ). 

The following proposition represents the essence of dynamics (see also [21 [8HT0] ) (in what 
follows, for the sake of brevity, we will not consider thermodynamics which, along with the 
motion equation, defines more fully the concept of mechanical system [H ITT]). 
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Fundamental principle of dynamics. For a mechanical tube Ac A there exist a Carte- 
sian frame So and a parameterization T of the r-image such that the vector fields f x and v x 
are solutions of the following 

j t V°(A t )=F°(A t ), A t C A, t€T (13) 

In this case 

1. the frame So and the parameterization T are called inertial (the frame is also called 
that of reference ); 

2. the aggregate a = {cr 3 , a t , yU, Vi 6 T, A f C A t , V(A t ), J-"(A t )} is called mechanical 
system; 

3. the set A t is called (actual) shape undergone by the mechanical system at t G T; 

4- the differentiable map T — > {A t , t G T} is called motion of the mechanical system /§[/; 

5. relation |73|) is called motion equation; 

6. the integral K.(A t ) is called scalar measure of motion of the mechanical system; 

7. the screw V(A t ) is called vector measure of motion of the mechanical system; 

8. the screw J^iAf) is called screw measure of impressed action of the mechanical system 
a e = {(73, a t , n, V£ 6 T, A' C A t , V(Af), — T(At)} on the mechanical system a (it 
defines the action of the environment A e on A). 



Relation (13) can be transformed in the vector form 

^P{A t ) = J-(A t ) (14) 

Note that the parameterization and the frame So are elements of Galilean group of trans- 
formations of A 4 which preserve intervals of time and the distance between simultaneous 
events [6]. 

The set A£ C At is called set of concentration of the measure M. if Ai(A t ) = for any set 



At C At \ A£. We shall assume that relation (14) is true only for A t C A£, t G T 



1.3. Concept of body 

The concept of a body is the subject of various mathematical formalizations. For example, 
one may represent a body as a point-wise set, an element of Boolean algebra, a differentiable 



manifold, a topological or measure space [TJ [5j [TT] where a map into the space of shapes is 
considered. But there is a small obstacle: we must also transfer masses and forces to body 
shares. 

If we do it in some way then the mathematical abstraction - body with mass and force - loses 
the primitive nature. In order to work out a mathematical theory of mechanics we have all 
the necessary: shares kinematical and dynamic structures attributed by them. That is why 
(if it is necessary) we shall use the following conventions in the case of a given mechanical 
system a = {a 3 , a u /i, Vt G T, A t C A t , V(A t ), ^(A*)}: 

1. the body is that takes some shapes A t C A t in 3-dimensional affine space at some 
instants of time (cf. Aristotel, Physics, III, 5, 204b); 

2. the change of shape undergone by a body from one instant to another is called the 
motion of that body (due to the principle of determinacy); 

3. the positive number A4(A t ) is the body mass; 

4- the screw J- (At) is the force impressed at the body. 

Figuratively speaking, it is just what we see on the movie screen, achieved by a range of 
consecutive film shots [11J. 

Note that according to Glossary, Earth Observatory, NASA: force is any external agent that 
causes a change in the motion of a free body, or that causes stress in a fixed body. 

While the concept of a mechanical system has strict mathematical sense, the concept of a 
body given above has only descriptive character, being a tribute of very fruitful tradition. 

1.4. Generalization of mechanical system concept 

The non-trivial nature of the mechanical system concept can be seen from the fact that we 



may postulate the following equation of motion (instead of (14)) 

d 



where 



d V(A t ) = F(A t ) + fi(A t ) + J- c (A t ) (15) 



V(A t ) = f I X J^^(dx) (16) 



is the inhomogeneous screw measure of motion of the mechanical system (it is not necessary 



to think that p = p x v x in (16)); J 7 is the inhomogeneous screw measure; the inhomogeneous 
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screw measure Ti is so called increment velocity of the measure V; the inhomogeneous screw 
measure J- c is so called constraint action. 

Assume that T c = J-'mt + J- e xt where Ti n t is formed by internal constraints of the set A t 
while J- ext is formed by external constraints. 

Note that the base point in mechanics is that the motion of bodies is caused by interaction 
with their environment in the universe A which is described in the terms of force, moment 
of force and 'innate' moment [T2] . However one must pay attention to constraints on bodies 
and their parts as well as interchange of masses, linear and angular momentums. 

1.5. Derivatives of some measures 

Let us present the set A t as the union of the set Af of the pure points entering into it, with 
the set A" c of its continuous points. We will assume that the last set has the surface <9A" C 
which is Lyapunov's simple closed one [T2] . 

Due to Gauss-Ostrogradsky (divergence) theorem, we have (see also [TT] ) 
d f d d 

-^M(At) = J X K «Aj t Px + pxdivV x )fiac(dx) + Y^i^PxJPppiXk) 

We shall assume that the function p x is defined by the continuity equation for continuous 
points §iPx + divp^ = v x and for pure points - j[ t p Xh = v Xk where u x depicts the generation 
(negative in the case of removal) per unit volume and unit time of the measure Ai. Terms 
that generate (y x > ) or remove [y x < 0) are referred to as 'sources' and 'sinks' respectively. 

In what follows, we shall assume that all vector and tensor sliders are homogeneous. 

According to [TTJ 



±V(At) = I X Ar (^ Px + l p ^v x )fi ac (dx) + J2(jJ Pxk )»M ( 17 ) 

* k 



As for continuous points (see also [TT] ) 



|-P* + P*div^ = P ^p* + ^r* 

at at 



from (15) and (17) follows 

/ X ^ Px Jt lVX + = ^ A *) + ^( A *) + ^( A *) ( 18 ) 
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2. Specifying mechanical systems 

Show how the given above axiomatics relates to the conventional mechanics. In the first 
place exemplify the notion of skew screw bi-measure. In the conventional mechanics it is 
considered that there is the gravitational interaction between bodies. It can be formalized 
in the following way. Let a skew screw bi-measure &(A t , A|) be such that 

*(A*, A t e ) = / X Aj 9 *pMdx), g x = 7 / xA*-y) PyKdy ) 



(x - y)\\ 3 

where 7 is a positive (gravitational) constant, the /z-integrable homogeneous slider p x l 9x is 
defined at x G A t . 

Then the screw Q(A t ) = f \P(A t ,Af) can be called measure of gravitating action of a e upon 
a [SJ. One may take this screw as the screw measure J~{At) of impressed action. 

Assume that the increment velocity of V is given by the following Lebesgue-Stieltjes integral 

Fi{A t ) = J x A / x p(dx), A t c A* (19) 

where fa is its density. 

Let no external constraint be. 

2.1. A mass— point 

Consider a world-line Ac A whose image in A 3 x T generates the curve {x(t) G A t , t G T}. 
Assume that the points x(t) are pure, i.e., x(t) = Xk{t). Then the mechanical system 
a = {03, a t , n, Vt G T, x(t) G A t , p x , u x , f x , ^} is called mass-point. 



From relation (18) follows that 



Px-^v x + v x v x = f x + £ x (20) 



If v x = and £ x = 0, then equation (20) is known as second Newton's law where f x is the 



impressed force acting at the point x = xt G A t with the mass M. k = p x \i\.. 



If v x 7^ and £ x = v x u x where u x is the velocity of mass gain or loss, then equation (20) is 
known as that of Meshchersky [10] . 

Remark 1. A classical example of mass-points with constraints is the mechanical system 
khown as pendulum. 
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2.2. A rigid body 

The mechanical system a p = {03, o~ t , Vt G T, A t C A t , Vx G A t , p x , / x , is called 
rigid body if 

1. the sets A t are bounded and closed; 

2. the constraints applied on its points keep distances between them not changing with 
time; 

3. the constraints are ideal f73] /. 

A rigid body may contain continuous and pure points [TO] . 

2.2.1. Newton— Euler equation. At any time instant t* consider the set A t *. Let a 
Cartesian frame E p be attached to the set under consideration. It is plain that the frame 
takes the same position in all sets A t . In the frame these sets are immobile, coincide one 
with another and form the set noted as A p in the frame S p . We shall say that the frame S p 
is attached to the rigid body a p . 

Let us bound the vectors v$ )P and uo tP at the point O p and the vectors v x and ljq )P at points 
x G E 3 . Then due to (|3j we have the field V 0tP = l u °,p> v °>p = {uj 0jP ,v , p + r XiP xu Qp , Vx G E 3 } 
known as kinematic slider. The coordinate representation of its reduction (twist) V^' 
co1{vq p , 0Jq p } G R.6 is called vector of quasi-velocities. 

Lemma 1. There is the following relation [5] 



rtw,p 



]V x ,wr,p fi)X-irtw,p p.x 

l p — W p V Qp , W p 



T _ r px 
p.x 

r px _( r px\2 
p,x v p,x) 



Proof. The statement is true as 



jv x ,wr,p 



I 


V P = 

u x 


I 












. p ' x . 



( V 0,p r p,x U Q,p) 



T _ r px 

p,x 

r PX —( r PX) 2 

p,x V p,xl 



0,p 



0,P, 



According to the rigid body definition the internal constraints are considered as ideal and 
thus [lj 

•Fc(A t ) = 



From relations ( 18 )— ( 19 ) follows that [TT 

d 
~dt 



" x J± ' U X l Q — l Q > l 



(21) 
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Hence we have 

/ x Ap \p^lj t (Lo, P 6 x p v^ p p ) + v x e;vZ' p Wx) = I x Ap [it* + it p Mdx) 

As the twist V Q ^' P does not depend on points x, the following statement is true. 

Theorem 2. The motion of a p (w.r.t. Sq in the frame S p ) is described by the (Newton- 
Euler j equation [5 J 

QXT + (Q P + *Z® P )K' P = ^ r,p (A P ) + JT P r ' p (K) (22) 

where Q p =Jx Ap @pV x K dx )> °p = ! X Ap O p Px^{dx). 



2.2.2. Systems of consecutively connected rigid bodies |14] . Let us consider a 
system of k + 1 consecutively connected rigid bodies a p , p = 0, k (the rigid body a is 
immobile). Its motion is depicted by the following Newton-Euler equation 

AV a + BV a = F a (23) 

where A and B are known matrices, V a = co1{Iq*p' p }, F a = co\{J r p r,p + J^ r ' p }, p — 1, k. 



Newton-Euler equation (23) is considered w.r.t. 'absolute' quasi-velocities Vq™' p of the rigid 



bodies (calculated in S p w.r.t. the main frame So). But in practice there are only the 
'relative' quasi-velocities V™f of the frame S p w.r.t. S v _\. Thus we must connect the 
'absolute' quasi-velocities with 'relative' ones. 

Lemma 2. For a system of consecutively connected rigid bodies there is the following com- 
position rule [5] 

v o7 = 2 L Z v s-C K-X = ( ^" 1,s ] (24) 



s=l 




where L% ts given as in {§) for L™ = C® S D^ S = (D^CfJ' 1 = L™f. 



Proof. As the rigid bodies are connected consecutively there is the relation Co jP = Co, s C StP . 
With differentiating it we have uj p p = ui p s + = co p s + C PtS cu^ p . 



Let us define the vectors d S)P = O s O p and do, s = OqO s , then d , P = OqO p = do, s + d StP , 

„,P — „,P _i_ Ap. ap —n as jp. _ q .p < fj- _ p _,,pxap — n,P _i_//px,.iP Hpnrp 

u 0,p ~~ u 0,s ' a s,pi a s,p — ^"s,;)) a s,p ~ u s,p * ^p,s a s,p ~ u s,p UJ s,p a s,p — u s,p ' u s,p UJ s,p- J-J-cllbB 

v o, P = v o,s + C P,sd s s * p u S s,p + v%, and 

TrtW,p TrtW,p , jtWTrtW,S 

V 0,p — V 0,s "T lj p,s V s,p 
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Hence we have (24). 

From ( 24 ) follows the equation of kinematics 



Va = LV r (25) 

where V a =co\{ V^' 1 , . . . ,Vq W p ' p , . . . ,V*J' k }, V r =co\{ V^'\ . . . ,V^ p , . . . ,V^ k }, L is the tri- 
angular matrix with blocks being functions of 'relative' frame rotations and translations 
(and their velocities). 

Thus we have 

ALV; + (AL~ +B)V r = F a (26) 
where U is analytically calculated due to relation u7\. 



It is easy to see that the matrices of relation (26) depend on rotation matrices (and linear 



and angular quasi-velocities, too) that is why equation (26) must be considered along with 
the Euler kinematical relation 




Fig. 1. Multibody system graphs 

2.2.3. Multibody systems with tree— like structure. Consider a multibody system 
with tree-like structure given by the graph in Fig. 1A. Let vertices j % represent the system 
bodies or the origins of the attached Cartesian frames £j where the index i numbers the tree- 
tops, the index j numbers the rigid bodies from the base to the corresponding tree-tops. In- 
troduce Vj^' 1 as quasi-velocities characterizing rotation and translation of the frames £j w.r.t. 
£™ . Then we have the sets { Vj 1 , V$ , V$ , V$ , V$ } , { Vjf , V$ , V$ , , , v£$ } , { Vj 1 , 
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yl,i yl,L yll yi,Z yZ,6, ryU,lyL,l yl,L ylfl, r yVl yl ,1 yi fi y^O | |yU,l yl,l yi,0 ^,61 

^1,2 > ,/ 2,3 J v 3,4 > "4,5 J v 5,6 J) \ v 0,l v l,2 > v 2,3 ) "3,4 J) I "0,1 "l,2 "2,3 "3,4 h I "0,1 J "l,2 ' "2,3 J "3,4 J 
j fT/0,1 t/1,1 t/1,1 t/1,1 t/1,1 -| TT/0,1 t/1,1 t/1,1 t/1,1 t/1,2 t/2,21 TT/O. 1 T /1,1 t /1,1 t/1,1 t/1,2 

^e?}, ^2, V?/, ^4 4 }> ^2 , Vof}, {Vj 1 , V^, V^i, 1/ 5 f} with the same sub- 

scripts as in the case of consecutively connected rigid bodies for the relative and absolute 



quasi-velocities. This case is considered above that is why we arrive at relation (25) with the 



known matrix L and V a = col{V °i\ Vq£, Vq£, V$, Vq£, V$, V^, Vq£, V '£, V^, V$, Vq£} 



rl.l 



'1,1 



1.1 



rl.l 



'1,2 



r2,2 Tr2,3 



-1,4 



'1,5 t/5,5 



r5,6 n 



t / lfT^U,! T 7"1,1 T/1,1 t/1,1 T/1,1 T rl,z T t2,£ T/2,3 T/1,4 T/l,t> T ro : o 

Remark 2. I7ie results obtained can be immediately applied to systems with loops, e.g., if 
in the system under consideration (see Fig. IB) the vertex 6 2 is connected with 6 3 by the 



edge (6 , 6 ). In this case relation (22) is the same, but in the case where constraints are 
considered there are the following additional constraints (5 2 , 6 2 ) + (6 2 , 6 3 ) + (6 2 , 5 2 ) = and 



✓^2,2 ^^2,3 y-y2,2 
L/ 5,6 ( - y 6,6 ( - y 6,5 



I. 



2.2.4. Parameterization of rotation matrices. The order of system (26)-(27) may be 



Ci 



,c 2 



,c 3 



(28) 



reduced. To this end one uses different parameterizations of rotation matrices. 
2.2.4.1. Euler angles. Let C p _i iP = CiC 2 C 3 where 

cos $ sin i? 

1 

- sin $ cos <£> 

are so called the simplest rotation matrices; ip, and ip are Euler angles [To] . 

Introduce the triple A p _i jP = col{<£>, #, ip} as a parameter. Then there is the matrix D p _i tP 
such that [5] 



1 

cos (p — sin 9? 

Sin If COS <£> 



cos i/j — sin ^ 
sin ip cos -0 
1 



(29) 



Hence equation (26) must be considered along with the following relation 



V — i p 

a p-i,p ~ ^p-i^p-x 

and C p _i )P = Cp_i jP (Ap_i )P ) if the matrix D~i iP exists. 



(30) 



2.2.4.2. Fedorov vector— parameter. To parameterize rotation matrices we may intro- 
duce Fedorov vector-parameter [T6] . 
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Definition 3. [16J The number triple f G E3 is called Fedorov vector-parameter of a 
rotation matrix C, if it corresponds to the following matrix 

r = (c - i)(c + 

The inverse map of Cayley restores the rotation matrix 

c=(i + d(/ - rr 1 

It is easy to be verified (for example, by means of Maple®) that the following relations are 
true 

r _ C-C T (i-||/f)z + 2// r + 2/* 

1 1 + trC" " 

Let the rotation matrices C p ^ k have Fedorov vector-parameter f Pjk . It is known that it is an 
eigenvector of C Pjk , i.e., C p , k f Pjk = f Pjk G E 3 . 

As the space E 3 has 3 bases e , e p and we may define the following vectors 

g P , k = E /a*. = E = E 

i i i 

where col{/i, / 2 , / 3 } = 

Definition 4. T/ie vector g Pik is called vector of Gibbs, while f Pjk is called vector of Ro- 
drigues. 

Remark 3. In fTbJ it is explicitly pointed out that Fedorov vector-parameter f P)k is Gibbs 
vector g p , k , and f Ptk = g Pjk as there is no other basis except the canonical one in [16]. In 
general the vector g Pyk does not coincide with f Pjk (see Fig. 2), as the bases eo, e p and e k 
are different. Moreover the name of vector is used here conditionally as there is not the 
parallelogram rule for the vectors of the kind g Pjk and r P)k (see also YT2\ [Tdj \PT^). 

It is easy to see that the vector f Ptk is collinear with the instantaneous angular velocity uj P)k , 
and thus it defines the rotation axis (uJ Pjk is the half of the vector of finite rotation [1TJ ) . 

There is the following relation [TT] 

2 

Wp-i, p — || j pvp-l,P + r P-l.P X r p-l,P> 
' II P— IjP II 

where x means vector product. 
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Fig. 2. Gibbs and Rodrigues vectors. 



As the vectors ujp-x tP and r p _i iP are collinear, we have relation (29) where A p _i iP = f p -i tP 

2 



and -Dp_i, p - 1+ y p _ ip p(I + fp-i, p )- 



Thus equation (26) must be considered along with the following relation 

4-1, - 1 a + u wx* - /;-i lP ) 2 (^ + (si) 

and C p _i )P = Cp_i )P (/p_i jP ). 

2.2.4.3. Euler— Rodrigues parameters. To parameterize rotation matrices we may use 
quaternions. 

Definition 5. TTie set A = {A G R, A G V 3 } is called quaternion. 
Quaternions generate the algebra with the quaternion product 

A o M = { A /io - (A, P) , Ao/t + /ioA + A x p} 

where M = {fj, ,p}. 

Any vector A can be imaged as a quaternion A with the zero scalar part. That is why we 
may define the quaternion product of two vectors A and p as follows 

A o ft — {— (A, p), A x p) 

There exists the unit quaternion A p _i iP = {Ao, A p _i iP } (with ||A p _i iP || = 1) such that [T8] 

A p _i, p oo; p _ 1)P oA p _i iP = cJ p _i )P , oJ p _i iP = — 2A p _i iP oA p _ lp , A p -i, p — ^ ^p-i,p°^v-x,p (32) 

where A p _ l p = {A , — A p _i iP } is conjugation of A p _ 1]P . 
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Cp-i,p(Ap_ lp ) 



(33) 



Let us denote col{wi, u 2 , w 3 } = f w£_ 1|P , col{Ai, A 2 , A 3 } = f AJ_ 1(P and col{A , Aj_ ljP } d = AJ_ 1|P 
then the orthogonal matrix C p _i iP corresponding to a rotation by the unit quaternion A p _i iP 
is given in the following form [T8] 

Aq + A^ — A 2 — A| 2A1A2 — 2A0A3 2A1A3 + 2A0A2 
2A1A2 + 2A0A3 Aq — A^ + A| — A3 2A2A3 — 2A0A1 
2A1A3 — 2A0A2 2A2A3 + 2A0A1 Aq — A^ — A2 + A3 

The quadruple A p _ lp is known as that of Euler-Rodrigues parameters. 

From ([32]) follows [IE] 

— Ui —u 2 —uj 3 

LO\ 0J 3 —OJ2 
OJ2 — UJ3 LO\ 

oj 3 uji —uj\ 



A P. 



(34) 



Hence equation (26) must be considered along with relations (32)-(34). 



From (33) follows also that there is the matrix -D p -i, p = -D p _i >p (A p _ lp ) such that relation 



p 

p-i.p 



is true. 



2.2.5. Lagrange equation of II kind. Let A p _i )P be a triple of Euler angles or Fedorov 
vector-parameter. Introduce the following notions. 

Definition 6. 1. The vectors q p -x iP = col{d p _ l p , A p _i iP } and 9 p _i iP = col{d p l liP , A p _ ljP } 
are called canonical generalized coordinates and velocities of the frame S p in the motion 
w.r.t. the frame £ p _i/ 
2. the relation 



K-tv = M p-hp1' P -i^ M p-i,p = diag{/, £> p -i,p} 
is called equation of kinematics of £ p -frame w.r.t. £ p -\. 



(35) 



From relations (26) and (29) follows the Lagrange equation of II kind 

A(q)q" + B(q,q')q- = T 



(36) 



where A{q) = I.'M'AL M. B(q,q-) = L T M T [ALM' + {AL + B)M], T = L T M T F a , M = 
diag{M p _i iP }, q = col{g p _i iP }. 

In the many cases there are constraints on motion of multibody systems, and the matrix iV 
exists such that the matrix N T N is non-degenerate and we may introduce the generalized 
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coordinate q = Nq G R m where the natural number m is not more 6k [19]. Then from 



relation (36) follows 

Aq" + Bq- = f 

where A, B and J 7 are known matrices and column. 

As to the quadruple A^_ x , we may replace A p _i ;P with Ap_ lp in the above definition and 



equation (36). It is clear that the corresponding matrix A proves to be singular. Under some 
assumption this equation is equivalent to a system of differential equations in Cauchy form 
and algebraic ones. The algebraic equations can be treated as constraints on the multibody 
system motion. It means that we may introduce 'new' coordinates, e.g., Euler angles or 
Fedorov vector-parameter, in order to obtain the Lagrange equation with a non-singular 
symmetric matrix A. 

2.3. A continuum 

2.3.1. Notion of continuum. We shall assume that the sets A t C A t are bounded and 
closed, all their points are continuous, and their surfaces dA t are Lyapunov's simple closed 
surfaces [12]. Let fi2 be the restriction of fi on the surface dA t , n x be the normal to this 
surface. 



Due to [20J constraints being in a small vicinity of x G At cause stress. Define the internal 
constraint action as follows JTJ [11] 

^mt(At) = Jxa A J Txnx ^{dx) 

where T x is called stress tensor. 

Due to Gauss-Ostrogradsky (divergence) theorem, we have [TT] 

Fint{A t ) = J x At div/ 7; / u(rfx) 

Take a point y(t) in a small vicinity of x(t) G A t at an instant t G T and define their 
radius-vectors f x (t) and r y (t) (in Sq) and the vector h(t) — r y — f x (t). Then there is the 
Cauchy-Helmholtz relation [11] 

v y (t) = v x (t) + ]r[dv x /dr x + {dv x /df x ) T } h(t) + ^[dv x /dr x - (dv x /dr x ) T ] h(t) 

where ^[dv x /df x + (dv x /dr x ) T ] is known as tensor of strain velocities; \[dv x / df x — (dv x / dr x ) T ] 
is known as spin-tensor at the point x G A t at the instant t. 
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Define the tensor S x (t) as the solution of the following equation 

$'x(t) = ^[dv x /df x + (dv x /df x ) T ] 
with initial data S x = I, t = to, X is the identity (spherical) tensor. 
The tensor S x is called strain one [5] . Let us define U x as S x or S' x . 

Definition 7. The mechanical system a = {a 3 , at, t 1 , Vi 6 T, A t C A t , Va; G A t , p x , v x , 
fx, £x> Tx] is called continuous medium or continuum of Navier-Stocks-Lame class if the 
tensor T x is an isotropic map oflA x , i.e., invariant w.r.t. orthogonal transformations. 

2.3.2. Quasi— linear isotropic matrix— functions. 

2.3.2.1. 3— dimensional case. For any 3 x 3-matrix U the aggregate PUQ is an isotropic 
function of U if the matrices P and Q are proportional to / with scalar coefficients being 
invariant w.r.t. rotations. 

Define the matrices 

E 1 = (tr£/)J, E 2 =U, E 3 =U T (37) 
where I is the identity matrix. 
Consider the following linear combination 

T = nEi + r 2 E 2 + r 3 E 3 (38) 

where are invariant w.r.t. rotations (they can be functions of the time, invariants of U 
and so on). 

Theorem 3. All isotropic quasi-linear 3 x 3-matrix functions of entries of U are given by 



relation (38) [21] 



2x2 matrix, / 



2.3.2.2. 2— dimensional case. Let U be 2 x 2 matrix. It is easy to see that for 2 x 2 
matrices P and Q the aggregate PUQ is an isotropic map of U if P and Q are of the kind 
rl + rl where the scalar coefficients r and r are invariant w.r.t. rotations, / is the identity 

-1 

1 

Introduce the following matrices E x = (tr[/)I, E x = (trIU)I, E 2 = U, E 2 = IU, E 3 = U T , 
E 3 = U T I, £ 4 = IU T , E 5 = UI, E 6 = IUI, and E 7 = IU T T. It is easy to see that there are 
6 linearly independent matrices, e.g., Ei, E x , E 2 , E 2 , E 3 , and E 3 . 
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Thus there is the set of isotropic quasi-linear 2 x 2-matrix functions of entries of 2 x 2-matrix 
U (invariant w.r.t. SO(R, 2)) 

T = n (tr U)I + n (tiIU)T+ r 2 U + r 3 U T + r 2 IU + r 3 U T I (39) 

where fj and are parameters being invariant w.r.t. orthogonal transformations. 



2.3.3. Symmetry of stress tensor. From equation (18) follows (in the inertial frame So) 

mm d 

Px-^v x + v x v x = p x g x + £ x + div%, % = Tj (40) 
We may introduce the following constitutive relations with the help of symmetrizing relations 



(38) and (39): in the 3-dimensional case 

% = r l + ri (tiU x )X+ r 2 U x 

and in the 2-dimensional case 

% = r l + ri (trZ4)X + r 2 U x + r 3 (lU x - U X T) 

where are rheological coefficients (parameters being invariant w.r.t. orthogonal transfor- 
mations); the tensor X corresponds to the matrix /. 

2.3.4. Correct continua. A continuum is called correct if the corresponding constitutive 
relation is invertible [SJ. With the help of routine calculations we see the following statements 
to be true: 

1. In 3-dimensional case let (3r*i + r 2 )r 2 ^ 0. Then there exists the inverse map 

U x = n l + rii {tv%)I + n 2 % 

where 

r ri 1 

n ° = i — i ' n i = -r, — j v n 2 = — 

3ri + r 2 r 2 (3r 1 + r 2 ) r 2 

2. In 2-dimensional case let (2r x + r 2 )(r| + 4r^) ^ 0. Then there exists the inverse map 

U x = n l + rix {ti%)X+ n 2 % + n 3 (Z% - %X) 



where 



r -rir 2 + 2ri r 2 — r 3 



2ri + r 2 (2ri + r2)(r 2 + 4rg) r 2 +Ar 3 r 2 +Ar 3 
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2.3.5. Kinds of continua. If IA X = S x and ro = the continuum is called elastic material, 
if U x = S x and r > (called Pascal pressure) the continuum is called viscous fluid [22J. 

The continua given above coincide with the continua used in continuum mechanics in the 
following cases [U [22] 
— the Pascal pressure r is positive and r 1 = r 2 = r 3 = (ideal fluid); 

ro is non-negative and r\tiX+ r 2 ^ — > trTx ^ — rotrX (correct continua) (here 
X is used as 2- and 3-dimensional identity tensors, respectively) . 

2.4. Systems with inhomogeneous screw measures 

Show how the systems with inhomogeneous sliders can be realized in the conventional me- 
chanics. 



2.4.1. Multiphase systems. Equation (15) is realized for multiphase systems (equations 



(6.34) and (7.11), given in . can be written in the form of (15)). Here the stress tensor 



proves to be non-symmetrical [UJ, and the motion equations are six-dimensional. 

2.4.2. Elements of Eulerian mechanics. The mechanical sense of the vectors p x and q x 



in equation (15) may be clarified in the framework of Eulerian mechanics [12] (we are not 
going to discuss its meaning as a base of mechanics). 

A particular case of bodies is the well-known mass-point being the fundamental concept 
of theoretical mechanics. It is considered as the unique model of a natural things having 
infinitesimal sizes, but possessing masses. Is such model a universal one? To answer this 
question, it is necessary to address to physics. The modern physics draws the following 
picture of a material objects: it consists of molecules, atoms, protons, neutrons, electrons, 
neutrinos or from their aggregates which are called clusters. What of these objects leads 
to the concept of a mass-point of theoretical mechanics? Let us take, for example, an 
electron. Its sizes are extremely small, it possesses some mass, so as though, it may be 
modeled as a mass-point. But here that disturbs us. It has appeared that at decoding 
and interpretation of tracks of nuclear particles, including electrons after their collisions, 
it is necessary to consider spins of these particles, to be exact, their angular momentums. 
Angular momentum is connected with rotation of these particles. But by definition a mass- 
point cannot rotate. It means that even a such small object as an electron cannot be modeled 
as a mass-point. Let us take a larger object, for example, a cluster or crystallite of some 
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polycrystalline metal. Certainly, it is possible to model motion of its center of masses as 
motion of some mass-point having the same mass, as well as the mass of cluster or crystallite. 
But a cluster or crystallite can rotate round the center of masses. Thus a cluster cannot be 
modeled as a mass-point, too. 

That is why Eulerian mechanics supplies points of the sets A t , t G T, with translation 
velocities v x and angular ones u x G V 3 , as well as densities A x , B x and C x of generalized 



inertia tensors. Then the kinetic energy (11) is introduced by its positive defined density 
k x d — h{v x ,A x v x ) + (v x ,B x u x ) + ^(uj x ,C x uj x ). After that one defines the vectors p x d = 
J-/c x = A x v x + B x uj x and q y d = r* x p x + q x where q x =^k x = B T x v x + C x u x is the density of 
so called dynamical spin. 



With the help of p x and q x the slider is introduced, equation (15) is postulated 



Realization of equation (15) in Euler mechanics is motion equations of point-bodies and 



their systems, thin rods and so on 



Conclusion 

In mechanics there is mainly absent the understanding that motion of bodies and interaction 
between them can be described with the help of screws as it is considered as conventional 
that '. . .being very attractive representation of a system of forces and rigid body motions 
with the help motors and screws, nevertheless it has no essential practical value. . .' |23j and 
that the screw calculus is not adapted for the description of continuum motion [24] . 

At the same time screw calculus gives useful, convenient and necessary tools which permit 
us to postulate the fundamental principle of dynamics in the differential form (see [5] and 



auhtor's paper 'On Foundations of Newtonian Mechanics', arXiv:1012.3633). However this 
form leaves in a shade many important features of rational mechanics that can be understood 
only with using the (stronger) local (primitive) integral form of the conservation (change) 
law for the vector measuare of motion (the differential form is applicable only in the cases 
where the divergence theorem is true). 

In order to obtain the integral form, the new notions of homogeneous and inhomogeneous 
vector and tensor slider-functions and screw measures are used, and the main mechanics 
measures, the equation of motion and the concept of mechanical system are introduced. It 
is shown that mass-points, rigid bodies, continua, multiphase systems, and point-bodies are 
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realizations of mechanical systems of the given axiomatics (see also [21 (TTJ [T2]). 
In this way we solve also the following problem: 

'. . . the dynamics of a continuous system must clearly include as a limiting case (correspond- 
ing to a medium of density everywhere zero except in one very small region) the mechanics 
of a single material particle. This at once shows that it is absolutely necessary that the pos- 
tulates introduced for the mechanics of a continuous system should be brought into harmony 
with the modifications accepted above in the mechanics of the material particle' |25j . 

The author would be highly grateful with whoever would bring any element likely to be 
able to make progress the development, and thus the comprehension, of the paper. Any 
comments, critiques, or objections are kindly invited to be sent to the author by e-mail. 
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